Magnitude and crystalline anisotropy of hole magnetization in (Ga,Mn)As 
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Theory of hole magnetization M c in zinc-blende diluted ferromagnetic semiconductors is developed 
relaxing the spherical approximation of earlier approaches. The theory is employed to determine M c 
for (Ga,Mn)As over a wide range of hole concentrations and a number of crystallographic orientations 
of Mn magnetization. It is found that anisotropy of M c is practically negligible but the obtained 
magnitude of M c is significantly greater than that determined in the spherical approximation. Its 
sign and value compares favorably with the results of available magnetization measurements and 
ferromagnetic resonance studies. 
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I. INTRODUCTION 

Following the discovery of ferromagnetism in zinc- 
blende Mn-based semiconducting compounds, 1 a theory 
has been developed that describes correctly a number 
of properties of those materials. 2,3 The theory, which is 
based on the Zener model of ferromagnetism mediated by 
band carriers, takes into account the complex structure 
of the valence band of zinc-blende semiconductors via 
the k ■ p method, and employs the molecular-field and 
virtual-crystal approximations to include the exchange 
interaction between the Mn d electrons and the valence 
band holes. 2,4 In particular, the magnetic moment of the 
hole carriers per unit volume, M c (p), has been calculated 
for a range of hole concentrations p. 2 The results, 2 cor- 
roborated recently by an independent calculation, 5 have 
shown that while a part of this magnetic moment can 
be attributed to spin polarization of the hole Fermi liq- 
uid, the orbital magnetic moment of the holes has also 
to be taken into account. However, in order to reduce a 
differential multi-band Landau equation to an algebraic 
one, it has been assumed in Refs. 2 and 5 that two rele- 
vant Luttinger parameters have the same value, 72 = 73- 
A quantitative error introduced but this approximation 
has been hard to estimate and, moreover, within such an 
approximation a significant warping of the valence band 
has been neglected. 

Motivated by the recent accurate magnetization mea- 
surements carried out by Sawicki and co-workers 5 and 
ferromagnetic resonance studies of Liu et al., 6 both com- 
pleted for high quality (Ga,Mn)As, we have decided to 
develop theory of hole magnetization in carrier-controlled 
ferromagnetic semiconductors valid for the arbitrary val- 
ues of 72 and 73. Our results show that in the experimen- 
tally relevant range of hole concentrations the theoretical 
values of M c (p) are much greater that those evaluated 
previously, 2,5 which allows us for a better description of 
the experimental findings. 5,6 Furthermore, we find that 
the crystalline anisotropy of M c (p) is practically negli- 



gible, so that the dipole interaction between the subsys- 
tems of Mn and hole magnetic moments does not con- 
tribute to magnetic anisotropy of (Ga,Mn)As. 

The calculation proceeds by the determination of the 
partition function of the holes occupying the Landau 
levels in in the presence of Mn spontaneous magnetiza- 
tion M, assumed to be parallel to the external magnetic 
field B. Hole magnetization is then obtained by differen- 
tiating the Gibbs thermodynamic potential with respect 
to B. In contrast to the case 72 = 73, a significant 
amount of symbolic calculations involving trigonomet- 
ric expressions is required, followed by an efficient nu- 
merical solution of a truncated infinite eigenvalue prob- 
lem. This methodology allows us to calculate and exam- 
ine hole magnetization for various crystallographic di- 
rections of the magnetic field and magnetization of the 
Mn-sublattice. 

The remaining part of this paper is organized as fol- 
lows: in Section II we present the details of the valence 
band structure of a zinc-blende semiconductor; in Sec- 
tion III we discuss quantized motion of an electron in 
uniform magnetic field; in Section IV we describe some 
technicalities of our algorithm; in Section V we present 
our results, and finally in Section VI we present conclu- 
sions of our work. 



II. THE k p MODEL OF THE VALENCE BAND 

Following classic works of Luttinger 7 as well as of Bell 
and Rogers 8 and Pidgeon and Brown, 9 the energy states 
in semiconductors of zinc-blende structure in a magnetic 
field (specifically InSb) have been considered in detail 
by Trebin, Rossler, and Ranvaud, 10 who obtained the 
effective mass Hamiltonian taking into account the two- 
fold degenerate conduction band of the symmetry Tg, the 
four-fold degenerate uppermost valence band of the sym- 
metry Tg, and the two- fold degenerate split-off valence 
band of the Ty symmetry. A model in which an upper 
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TABLE I: Matrices for the cross space of the valence and split-off valence band states (Ui = T 1 /, Uij = T^ ). 



Fg and conduction bands are explicitly included was 

developed by Pfeffer and Zawadzki. 11 

Here, we consider explicitly only the valence band, 

namely the bands of the and symmetry, denoted 

respectively as v and s, for which we choose the ba- 
sis C|3 3\ 13 1\ 13 _1\ 13 _3\ II 1\ II 14 
S1S V|2'2/'|2'2/'|2' 2/ 'I 2' 2/'l2'2/'l2' 111' 

Therefore, in our k ■ p model of the valence band of a 
zinc-blende semiconductor, the band electron Hamilto- 
nian has the block form 



H = 



(1) 



where 



W = ~-{\lik 2 ~ l2 [{Jl-\j*)kl + c.p.] 
m I 2 3 

- 27 3 [{ J x , J y }{k x , k y } + c.p.]} 



--k(J.B), 

m 



H ss = -(A 



2m 



lik 2 ) 



--(-+ T )(--S), 

n vs = —[-3 l2 (u xx kl + c. P .) 

- Ql3(U X y{k X} ky} + C.p. 

2 m 2 



(2) 



(3) 



(4) 



Here, {A,B} = \{AB + BA) and c.p. denotes cyclic per- 
mutations. The matrices Ji are the spin-| angular mo- 
mentum matrices, <ii are the spin-i angular momentum 
matrices multiplied by 2 (Pauli matrices in our basis), 
and Tj, are the cross space matrices introduced in 
Ref. 10, in our basis given in Table I (which corrects one 
misprinted sign in Ref. 10, Table I). Notice the Zeeman 
terms of the form given in Ref. 2 that properly accounts 
for the Lande factor of the free electron, go, a correction 
in the overall sign of H vs , as well as the factor of two 
in the first (kinetic) term of 7i vs that has changed with 
respect to Ref. 10. 

Now, to complete our model with the exchange interac- 
tion between the valence band carriers and the localized 
ci-electron spins, treated within the molecular-field and 
virtual-crystal approximations, we augment our Hamil- 
tonian with the p-d exchange matrix H p d, 



n 



pd 



(2(J-w) 6(U-w)\ () 



where the Bq parameter is proportional to the Mn mag- 
netization M and the p-d coupling constant /?, 



B G 



f3M 



(6) 



where Ap w 1.2 is the Fermi liquid Landau parameter de- 
scribing effects of the hole-hole exchange interaction, 2,3 
and w is the versor pointing in the direction of the 
Mn magnetization; see Ref. 2, Eqs. A13 to A18. 



III. QUANTIZED MOTION OF AN ELECTRON 
IN UNIFORM MAGNETIC FIELD 

The differential operators fcj in our k ■ p model are 



d eA a 
ox a n 

and satisfy k x k = eB/iU, or 

eB a 



£ a f3~fkf3 ky 



ih 



Following Luttinger 7 , by 



(fci - ik 2 ), 
(fci + ik 2 ), 



(7) 

(8) 

(9) 
(10) 



we introduce a pair of operators (a, a 1 *) satisfying the 
canonical commutation relation, [a, aJ] = 1. Here, s = 
eB/h, and the coordinate system "1", "2", "3" is such 
that B is along the "3" direction. Then, fc# = fc 3 (the 
momentum component along the direction of the mag- 
netic field) commutes with a and a) . We will express the 
effective mass Hamiltonian in terms of {a,a\ kn)- 

In the Landau gauge, the eigenvalues of n = a^a, to- 
gether with fci and fc 3 , number the energy levels (Landau 
levels) for a free electron in the uniform magnetic field B. 
The corresponding wavefunctions are given by 12 



ip = exp[i(fcixi + k 3 x 3 )]u n (x 2 - x 2 ), 



(11) 



where u n (y) = (s/tt) 1 ^ exp(— sy 2 /2)H n (y/sy)/\/2 n n\ is 
the wavefunction of a harmonic oscillator, and x 2 — fci / s. 
Since the vector potential A(x) = (—Bx 2 ,0,0) docs not 
depend on x\ and x 3 , we use periodic boundary condi- 
tions for those two variables. As far as x 2 is concerned, 
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the particle is localized by the magnetic field, therefore 
the boundary conditions can be ignored (at least in the 
thermodynamic limit, i.e. when the magnetic length is 
much smaller than the size of the system). On the other 
hand, the length of the system in the "2" direction deter- 
mines the range of integration over k\ . Hence, we have 
the formula (B3) of Ref. 2: 



G c = -k B T 



eB 
2irh 



dk 



H 



2tt 



2^ log{l + exp[ 1 — 7? ]}, 



i=0 



k B T 



(12) 



where we have changed the sign under exponent to reflect 
the fact that our particles are holes, while £i(fcj?) is the 
i-th energy level for an electron with k^ = ku (the re- 
sult differs by an additive constant corresponding to the 
valence band fully occupied by electrons). 



IV. DETAILS OF THE ALGORITHM 

Calculation of G c according to Eq. (12) requires precise 
knowledge of the energy levels In contrast to the 

case 72 = 73, the eigenproblem for energy levels does 
not decompose into 6x6 blocks. Therefore, we have 
to generate a truncated matrix of the Hamiltonian, and 
select from its eigenvalues only those that approximate 
eigenvalues of the infinite matrix. Our Hamiltonian is of 
the form H — hi(3l + h n (3n+h a (3a+hj tgia 1, + h a 2®a 2 + 



h, 



,t 



a* , where hi, h n , h a , and h a 2 are 6x6 complex 



matrices with entries being trigonometric expressions in 
the angles (6, <j>) specifying the direction of the magnetic 
field. 

Instead of just truncating the operators (a, a* 1 ) at 
some n max , we employ a procedure that is accurate 
for 72 = 73 in the sense that the generated matrix 
of the Hamiltonian is truncated at the non-zero 6x6- 
block boundary. This requires applying to H a uni- 
tary transformation H — > TZHTZ -1 , where 1Z = 1Z(6,4>) 
is a 6 x 6 unitary matrix that implements a rotation 
from the crystal back to the "1", "2", "3" coordi- 
nate system, so that in the spherical case the result- 
ing matrix does not depend on [0,4>). Then, since in 
the spherical case the wavefunctions assume the form 
(ciu n , c 2 u n+ i, c 3 u„ +2 , c 4 u„ +3 , c 5 u n+ i, c 6 u n+2 ), where c, 
are the unknown components of the eigenvectors, from 
the transformed 6(n max + 1) x 6(n max + 1) matrix we 
drop columns and rows numbered —18, —12, —11, —8, 
—6, —5, —4, —2, —1 (the sign minus means counting 
from the right/bottom), so as to preserve a number of 
whole 6x6 blocks in which (in the spherical case) the 
Hamiltonian matrix is non-zero. 

The generated matrix is diagonalized using the LA- 
PACK algorithm for band Hermitian matrices (only the 
eigenvalues are computed) 15 . An eigenvalue is then se- 
lected as correct if for two subsequent sizes of the trun- 
cated matrix we obtain eigenvalues that are the same 
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FIG. 1: Hole liquid contribution to volume magnetization, 
M c (p) as a function of the hole density p for 72 = 73 ( "spheri- 
cal approximation", Refs. 2, 5) and for the real band structure 
parameters of GaAs ("no spherical approximation"). The 
spin splitting corresponds to the saturation value of magneti- 
zation for Mn concentration x = 0.05. The calculations have 
been carried out for Mn magnetization along three principal 
crystallographic directions. 
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FIG. 2: The curves of Fig. 1 at low hole densities assuming 
various directions of Mn magnetization. 



within the numerical precision (this is called the "second 
truncation" test in Ref. 8). Great care must be under- 
taken because the set of selected eigenvalues does not 
have to be complete unless n max is large enough, i.e. 
there may be some lacking eigenvalues in the range be- 
tween the smallest and the largest eigenvalues selected. 
The minimal selected eigenvalue does not monotonically 
decrease with increasing n max in such a case. That is 
why we choose n max w 5000 as a starting point e.g. for 
the [110] direction of the magnetization M. 
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V. RESULTS 

We have performed computations for (Ga,Mn)As 
adopting the previously employed values of the band 
structure parameters and p-d exchange integral. 2 Fig- 
ure 1 presents the determined values of the hole mag- 
netization for three principal crystallographic directions 
of magnetization, compared to the data obtained for 
the case 72 = 73 = 2.58, which reproduce the ear- 
lier results. 2 ' 5 The value of the spin splitting parameter 
Bo = — 30meV has been assumed, which corresponds 
to an effective Mn concentration x e g = 0.05 if magne- 
tization is saturated. As seen, the curves for the three 
crystallographic directions [001], [110] and [111] of Mn 
magnetization overlap and are indistinguishable. This 
demonstrates that a strong warping of the valence band 
does not result in anisotropy of the hole magnetization 
in the experimentally important range of the hole con- 
centrations p. A magnified view of the region near p = 
is presented in Fig. 2 to show that in this range a small 
anisotropy of M c = becomes visible. 

The dependence of the carrier magnetization M c on 
the Mn magnetization M is presented in Fig. 3. Three 
curves for different values of the parameter Bq (—10, 
—20, and —30 mcV) arc displayed. For comparison, 
the data for Bq = lOmeV and 72 = 73 = 2.58 are 
also shown. An important aspect of the results de- 
picted in Figs. 1 and 3 is that the approximation 72 = 
73 = 2.58 employed previously underestimates signifi- 
cantly the magnitude of the hole magnetization. An ac- 
curate comparison between experimental and theoreti- 
cal results is, however, somewhat hampered by a rela- 
tively large uncertainty in the experimental value of M c , 
whose determination from magnetization measurements 
requires accurate information on the effective Mn con- 
centration. Nevertheless, as an example we consider two 
samples with x — 0.05, i.e. Bq — — 30meV, and hole 
concentrations p = 4.4 • 10 20 and 8.4 • 10 20 cm~ 3 . For 
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FIG. 3: Hole magnetization M c (p) for various magnitudes of 
valence band exchange splitting (Bg equal —10, —20, and —30 
meV). 



Sample No. 1 2 3 





Experiment 




n h -MF [cm- 3 ] 


1.24 ■ 10 20 1.48 • 10 20 


1.64 • 10 2() 




1.92 ±0.04 1.87 ±0.03 


1.80 ±0.02 




Theory 




V 


0.836 0.828 


0.822 


Mc/^b [cm -3 ] 


-1.46 ■ 10 20 -1.69 • 10 20 


-1.84 ■ 10 20 


9eS 


1.90 1.89 


1.88 



TABLE II: Comparison of effective Lande factor determined 
by ferromagnetic resonance (Ref. 6) to theoretical values ob- 
tained in the present work. 



these samples the previous model predicts M c = —0.25 
and —0.32 hb per one Mn ion, respectively. The cur- 
rent theory leads to M c = —0.36 and —0.49 \ib- These 
sets of results can be compared to the experimental data, 
M c = —0.36 and —0.80 /Ub, which we have obtained by 
taking a mean value of M c , as measured for x = 0.045 
and 0.056 before and after annealing. 5 We see that the 
present theory describes better the experimental findings 
but, as mentioned above, more quantitative comparisons 
requires detail information on the concentration of Mn 
spins that are not compensated by the antifcrromagnctic 
coupling to interstitial Mn neighbors. 

Another relevant experiment is ferromagnetic reso- 
nance, which leads to the magnitude of the effective 
Lande factor systematically smaller than 2, an effect 
taken as evidence for the hole contribution to the mag- 
netization dynamics. 6 A comparison of our calculations 
with the experimental data of Ref. 6 is given in Table II. 
The theoretical values of g e g have been calculated assum- 
ing riMn = 1-01 • 10 21 cm~ 3 , Bq = — 30meV, and using 
the formula similar to Eq. 2 of Ref. 6: 

Sg M nn M n + M c j[i B 

= 6n Mn ± sn h V, (13) 

9eS 

where S = 5/2, <?Mn = 2.0, s = 1/2, and V is the hole liq- 
uid spin polarization defined in Eq. 11 of Ref. 2. We see 
that our theory describe satisfactorily the deviation of the 
value of the effective Lande factor from 2, although for 
the sample with the highest Tq the deficit of the magne- 
tization is still larger than our calculation predicts. This 
disagreement may result from an underestimation of the 
experimental value of the hole concentration that was 
calculated based on the value of To 6 



VI. CONCLUSIONS 

In our work, presumably for the first time, an efficient 
approach has been proposed allowing for the calculation 
of the Landau level energies in the valence band within 
the k ■ p model of zinc-blende materials for an arbitrary 
direction of the magnetic field and non-zero hole momen- 
tum fcn taking effects of cubic anisotropy into account. 
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The model has been employed to evaluate the contri- 
bution of kinetic and spin energies to the volume mag- 
netization of the holes in (Ga,Mn)As. Surprisingly, the 
evaluated crystalline anisotropy of the carrier magneti- 
zation in the relevant range of the hole densities is neg- 
ligibly small. At the same time we have found that the 
magnitude of M c is about 50% larger than that calcu- 
lated earlier in the spherical approximation. We have 
demonstrated that our improved theory describes better 
the available results of magnetization measurements 5 and 
magnetic resonance 6 studies. 
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